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Abstract 


An  above-groiind  2.44  m  x  2.44  m  x  20  cm  rolled  homogeneous  armor  (RHA) 
plate  with  a  total  mass  of  9000  kg  located  asymmetrically  at  a  fixed  horizontal  61 
cm  offset  and  a  vertical  46  cm  stand-off  above  a  shallow-buried  mine  with  an 
equivalent  charge  weight  of  9.05  kg  TNT  and  15  cm  depth  of  overburden  in  dry 
and  wet  sand  has  been  modeled  in  3-D  using  the  CTH  hydrodynamic  code.  The 
purpose  is  to  compare  the  loading  and  response  of  the  anvil  plate  due  to 
detonation  of  a  mine  buried  in  dry  and  wet  sand  using  a  3-D  asymmetric  model 
of  a  mine-soil-hull  floor  interaction  problem.  The  model  allows  parametric  study 
of  complex  explosive-soil-structure  interaction  effects  due  to  variation  of 
explosive  types  and  plate  dimensions,  standoff,  types  of  soil  and  height  of  soil 
overburden  and  computes  momentum  loading  on  the  ballistic  pendulum  anvil 
plate  due  to  various  off-axis  locations  of  the  mine  in  dry  and  wet  soils  relative  to 
the  hull  floor  plate  located  on  the  underside  of  an  armored  vehicle. 
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ABSTRACT 


An  above-ground  2.44  m  x  2.44  m  x  20  cm  rolled  homogeneous  armor  (RHA) 
plate  with  a  total  mass  of 9000  kg  located  asymmetrically  at  a  fixed  horizontal 
61  cm  offset  and  a  vertical  46  cm  stand-off  above  a  shallow-buried  mine  with 
an  equivalent  charge  weight  of  9.05  kg  TNT  and  15  cm  depth  of  overburden  in 
dry  and  wet  sand  has  been  modeled  in  3-D  using  the  CTH  hydrodynamic 
code.  The  purpose  is  to  compare  the  loading  and  response  of  the  anvil  plate 
due  to  detonation  of  a  mine  buried  in  dry  and  wet  sand  using  a  3-D  asymmetric 
model  of  a  mine-soil-hull  floor  interaction  problem.  The  model  allows 
parametric  study  of  complex  explosive-soil-structure  interaction  effects  due  to 
variation  of  explosive  types  and  plate  dimensions,  standoff,  types  of  soil  and 
height  of  soil  overburden  and  computes  momentum  loading  on  the  ballistic 
pendulum  anvil  plate  due  to  various  off-axis  locations  of  the  mine  in  dry  and 
wet  soils  relative  to  the  hull  floor  plate  located  on  the  underside  of  an  armored 
vehicle. 


INTRODUCTION 

Both  armored  persoimel  carriers  and  light  combat  vehicles  are  increasingly  being  used  in  a 
support  role  for  other,  more  heavily  armored  combat  vehicles,  all  of  which  are  being  subjected  to 
much  greater  risk  from  a  variety  of  highly  lethal  antitank  land  mines.  As  a  result,  there  is  a  need  for 
modeling  and  understanding  the  interaction  of  mine  blast  products  with  structures  and  the  resulting 
loading  and  damage  mechanisms  inflicted  by  explosive  blast  and  impact.  This  understanding  is 
required  both  for  damage  assessment  and  protective  hardening  of  both  wheeled  and  tracked  vehicles. 
This  has  been  a  topic  of  considerable  interest  to  the  U.S.  Army  Research  Laboratoiy  (ARL)  at 
Aberdeen  Proving  Ground,  MD. 

An  interesting  study  [  1  ]  involved  collation  of  data  from  previous  scaled  experiments  related  to 
the  encounter  of  buried  land  mines  with  various  targets  including  flat  plates.  A  correlation  function, 
based  on  dimensional  analysis,  for  the  total  impulse  delivered  by  the  land  mine  and  overburden  on  a 
target  was  presented.  Vulnerability  of  tank  bottom  hull  floorplates  subjected  to  shallow-bmied  mine 
blast  was  studied  by  Norman  [2]  and  Hoskins  et  al.  [3] ,  who  employed  an  approximate  energy  method 
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approach  and  equated  the  strain  energy  absorbed  by  the  plate  to  the  energy  delivered  by  the  blast  to 
the  armor  plate  based  on  the  law  of  conservation  of  energy.  Various  simple  models  of  explosively 
driven  metal  were  developed  by  Dehn  [4]  and  compared  with  experiments. 

In  another  investigation  [5],  the  DORF  hydrocode  was  used  to  generate  the  loading  history  on 
a  target  plate  subjected  to  a  blast  load  from  a  land  mine.  The  mine  blast  loading  was  used  in  a 
structural  response  code  to  compute  the  response  of  a  clamped-edge,  square  rolled  homogeneous 
armor  (RHA)  plate.  The  predicted  deflection  and  total  impulse  were  in  satisfactory  agreement  with 
values  computed  using  an  empirical  correlation  function  derived  from  field  experimental  data. 

Dynamic  response  analysis  of  a  rectangular  simply  supported  plate  subjected  to  an  explosive 
blast  was  analyzed  by  Gupta  [6]  for  the  linear  elastic  case  and  compared  to  a  closed  form  solution. 
The  analysis  was  extended  to  the  elastic-plastic  domain  by  Gupta  et  al.  [7],  who  employed  the 
isotropic  bilinear  hardening  material  model  available  in  the  ADINA  nonlinear  finite  element  analysis 
code.  This  study  was  further  extended  by  Gupta  et  al.  [8]  to  the  response  of  a  generic  vehicle  floor 
model  subjected  to  empirically  obtained,  triangular  overpressure  loads  from  the  detonation  of  a 
shallow-buried  mine.  Although  qualitative  agreement  with  experiment  for  peak  deflection  could  be 
obtained,  the  residual  plastic  deformation  could  not  be  compared  effectively  due  to  lack  of  soil 
loading  arising  from  blast-soil-structure  interaction,  which  substantially  affected  the  estimated  blast 
loads  and  the  resulting  structural  response.  Another  recent  study  of  dynamic  elastic-plastic  stress  and 
deformation  response  of  a  generic  vehicle  floor  model  to  coupled  blast  and  impact  loads  from  a  plate 
mine  was  reported  by  Gupta  [9].  He  provided  an  inexpensive  method  of  evaluation  of  the  structural 
integrity  of  modem  vehicles  subjected  to  severe  spatially  varying  transient  loads  in  extreme 
environments. 

A  number  of  studies  have  been  performed  in  the  general  area  of  blast  response  to  stractures 
over  the  years.  However,  the  loading  mechanisms  from  blast-soil-structure  interaction,  such  as  those 
occurring  from  detonation  of  a  mine  buried  in  saturated  soil  resulting  in  moisturized  soil  are  poorly 
understood  at  present.  The  current  investigation  is  an  attempt  to  fill  this  void  and  to  simulate  loading 
of  the  anvil  plate  during  blast-soil-plate  interaction.  This  plate  acts  as  a  momentum  trap  for  the 
impulse  delivered  by  the  mine  detonation.  This  loading  is  needed  for  the  design  of  the  guide  shaft  for 
the  Vertical  Impulse  Measurement  Facility  (VIMF),  which  is  currently  under  construction  at  APG. 


OBJECTIVE 

The  primary  objective  of  this  investigation  is  to  simulate  detonation  of  a  mine  buried  in  sand 
bracketed  between  two  extremes  of  moisture  content: — dry  vs.  fully  saturated.  The  mine  detonation 
causes  propagation  of  the  blast-soil  wave  upward,  resulting  in  blast-soil-stracture  interaction  and 
consequent  loading  of  the  anvil  plate.  The  study  is  needed  to  obtain  improved  loading  functions  for 
combined  blast  and  soil  loading.  The  effect  of  soil  moisture  content  upon  the  momentiun  delivered  is 
one  of  the  important  questions  to  be  answered.  These  loading  data  can  then  be  transferred  to  a 
nonlinear  finite  element  analysis  code  to  obtain  transient  structural  response.  The  predicted  results 
can  then  be  compared  to  data  from  experiments  designed  to  capture  and  measure  the  loads  produced 
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in  these  events.  Such  capability  will,  in  turn,  enable  the  development  of  models  of  expansion  of 
blast/soil  events  in  which  moisture-dependent  soil  properties  play  a  significant  role. 

Accuracy  of  predictive  models  is  essential  to  the  development  of  innovative  vehicle  designs 
that  have  improved  survivability  characteristics  and  are,  at  the  same  time,  light  enough  to  ensure  rapid 
deployment  and  desired  load-carrying  capacity.  The  ultimate  objective  is  to  improve  survivability  of 
the  crew  in  vehicles  subject  to  land  mine  attacks. 


PROBLEM  CONFIGURATION 

The  problem  involves  the  worst-case  scenario  of  a  shallow-buried  cylindrical  mine  embedded 
all  around  in  two  different  soil  types:  dry  sand  and  wet  soil  with  a  15-cm  soil  overburden.  The  mine 
contains  explosive  with  an  equivalent  charge  weight  of  9.05  kg  TNT  and  is  asymmetrically  located 
with  a  horizontal  centerline  offset  distance  of  61  cm  relative  to  a  square  plate  244  cm  x  244  cm.  Plate 
thicknesses  of  1 5  cm  and  20  cm  were  modeled.  The  vertical  standoff  between  the  bottom  plate  surface 
and  the  top  of  the  soil  overburden  was  46  cm  as  shown  in  Figure  1.  The  cylindrical  configuration  of 
the  mine  has  a  diameter  of  32  cm  and  a  height  of  9  cm.  The  mine  is  located  inside  a  block  of  soil  366 
cm  in  diameter  with  a  total  height  of  124  cm. 


ASSUMPTIONS 

The  asymmetric  configuration  of  the  mine  as  defined  in  the  previous  section  is  a 
three-dimensional  (3-D)  problem  requiring  the  generation  of  a  3-D  model  that  may  be  expensive  to 
run  using  a  hydrodynamic  code.  However,  simplification  of  the  3-D  problem  to  a  2-D  model  using 
two-dimensional  (2-D)  plane  strain  assumptions  does  not  work  since  2-D  results  apply  to  a  semi¬ 
infinite  model  where  both  the  anvil  plate  and  the  explosive  geometry  are  extended  to  infinity  in  the 
Z-diection  contrary  to  the  finite  geometry  of  the  real  3-D  problem.  The  amount  of  explosive  and  the 
plate  mass  for  the  2-D  plane  strain  problem  are  significantly  higher  than  those  for  the  finite  3-D  model. 
Additionally,  loading  and  response  predictions  fi-om  CTH  computation  for  the  2-D  plane  strain  model 
cannot  be  easily  scaled  or  correlated  to  the  corresponding  results  fi'om  the  full  blown  3-D  analysis. 
As  a  result,  in  the  current  analysis  a  finite  3-D  model  was  employed  with  fixed  geometry  of  the  plate, 
the  soil  and  the  cylindrical  mine  close  to  actual  dimensions  to  be  employed  in  a  test  to  simulate  the 
worst  case  scenario. 

All  empty  spaces  between  the  plate  and  the  soil  block  and  all  around  these  configurations  were  filled 
with  dry  air  in  the  hydrocode  model.  The  wet  sand  is  assumed  to  be  fully  saturated  due  to  high 
moisture  content  completely  filling  the  void  volume,  which  is  approximately  27.4%  of  the  entire  soil 
volume.  This  ensures  a  worst-case  scenario  With  significant  contribution  to  the  resultant  momentum 
produced  in  the  momentum  trap.  In  the  second  case,  dry  sand  with  no  moisture  content  is  considered 
as  the  burial  medium  for  the  explosive.  This  allows  bracketing  of  the  influence  of  large  variation  of 
basic  soil  properties  due  to  maximum  possible  variation  of  moisture  content  upon  structural  loading 
and  response  for  the  most  commonly  available  soil  type. 
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Figure  1.  Initial  Configuration  of  the  Asymmetric  Mine-Soil-Structure  Interaction  Problem 

MODEL  DESCRIPTION 

A  3-D  asymmetric  model  with  a  61 -cm  centerline  offset  of  the  anvil  plate,  as  shown  in 
Figure  1 ,  was  generated  for  the  CTH  code  [  1 0, 1 1  ] .  A  3-D  Eulerian  mesh  was  overlaid  using  uniform 
relatively  coarse  grid  spacing  in  horizontal  X,  vertical  Y  and  lateral  Z  directions  throughout  the 
computational  space.  Each  cell  size  was  a  cubic  mesh  2  cm  x  2  cm  x  2  cm,  and  the  total  number  of 
cells  was  260,000.  An  aspect  ratio  of  1 .0  was  maintained  in  computational  space  to  ensure  numerical 
accuracy.  The  problem  region  was  extended  approximately  1 10  cm  vertically  above  the  top  plate 
surface  and  horizontally  beyond  the  left  and  right  geometric  edges  to  minimize  reflection  from  the 
computational  boundaries.  This  large  volume  of  space  filled  with  dry  air  accommodated  large 
displacements  and  tilt  angles  of  the  anvil  plate  expected  to  occur  due  to  asymmetric  loading  from 
mine  detonation.  Computed  prediction  of  peak  vertical  displacements  and  rotation  will  be  useful  in 
designing  vertical  guides,  guide  bars,  and  springs  for  the  experimental  test  fixture. 
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BOUNDARY  CONDITIONS 


The  steel  plate  is  assumed  to  be  freely  suspended  46  cm  above  the  top  surface  of  the  soil  block 
in  which  the  mine  is  embedded  and  is  initially  at  rest.  The  bottom  horizontal  soil  boundary  100  cm 
below  the  mine  is  assigned  a  transmitting  boundary  condition  where  mass  may  flow  in  and  out  of  mesh 
but  no  reflection  is  allowed  from  the  boundary. 

Both  left  and  right  vertical,  near  and  far  lateral  as  well  as  the  top  horizontal  boundaries  of  the 
computational  mesh  are  assigned  extrapolated  pressure  boundary  conditions  where  mass  is  not  allowed 
to  flow  into  the  mesh  but  the  pressures  in  cells  at  the  boundary  are  computed  using  extrapolated  values 
from  two  adjacent  cells.  This  boundary  condition  is  mostly  transmissive  with  some  reflection 
permitted. 


MATERIAL  MODELS 

The  model  involves  four  different  material  types.  The  first  package  in  the  problem  setup  represents 
the  TNT  explosive  for  the  land  mine.  The  next  package  represents  the  rolled  homogeneous  armor 
(RHA)  steel  plate.  This  is  followed  by  the  soil  package  surrounding  the  explosive.  The  final  package 
refers  to  dry  air,  which  is  used  to  fill  all  remaining  empty  space  between  and  around  the  previous  three 
packages  and  extends  up  to  the  edges  of  the  computational  mesh. 

The  default  SESAME  tabular  EOS  available  in  the  CTH  [11]  material  library  was  used  for 
RHA  steel.  This  EOS  treats  the  polymorphic  and  melting  ransitions  in  iron,  and  it  gives  good 
agreement  with  shock  wave  and  penetration  experiments.. 

Initially,  the  Johnson-Cook  [12]  constitutive  model  for  the  deviatoric  response  behavior  of 
steel  was  used,  which  resulted  in  appreciable  transverse  bending  deformation  of  the  steel  plate.  To 
reduce  large  deformation  response  of  the  steel  plate,  an  elastic-perfectly  plastic  constitutive 
relationship  with  a  yield  strength  of  1.5E09  Pa,  Poisson’s  ratio  of  .285  and  the  Johnson-Cook  fracture 
model  with  a  fracture  strength  of  -3.8E09  Pa  was  used.. 

The  material  model  employed  for  the  TNT  explosive  was  the  default  JWL  (Jones- W^kins-Lee) 
equation  of  state  available  in  the  material  library.  For  the  detonation  of  the  explosive,  a  programmed 
explosive  bum  routine  designated  as  HEBURN  in  the  hydrocode  was  employed.  The  initiation  was 
chosen  to  start  at  the  bottom  central  location  on  the  vertical  symmetry  axis  of  the  mine,  and  pressure 
contour  plots  were  generated  at  a  regular  interval  of  10-cm  radius.  The  detonation  velocity  was 
6.93E05  cm/s  [13].  The  default  SESAME  tabular  EOS  [1 1]  was  used  for  air  which  filled  all  empty 
space  in  the  computational  model. 

The  soil  material  models  were  rather  difficult  to  select  because  of  limited  information  available 
on  dynamic  soil  properties  of  widely  differing  soil  types  at  high  strain  r  ates.  Soil  properties  appear 
to  vary  considerably  depending  upon  the  location,  type,  moisture  content,  porosity,  aggregate 
size,  ambient  temperature,  and  humidity  conditions,  etc.  These  property  variations  can  significantly 
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affect  the  interacting  loading  mechanisms  and  generated  load  functions  as  well  as  the  resultant 
structural  response.  Most  dynamic  soil  equations  of  state  are  devoted  to  high  compression  regimes 
appropriate  in  underground  tests,  and  very  little  information  is  available  for  soil  models  in  tension  or 
dilatation  as  occurs  in  shallow-buried  mine  explosions.  Two  different  soil  models  i.e.,  completely 
dry  sand  and  fully  saturated  wet  sand  were  included  for  this  investigation. 

(a)  Dry  Sand 

Sand  is  a  porous  aggregate  of  particles  of  various  minerals,  most  commonly  silica  (fused 
quartz).  The  properties  of  sand  can  vary,  depending  on  the  degree  of  compaction,  the  types  of 
minerals  present,  and  the  moisture  content.  In  this  study,  fused  quartz  is  the  only  mineral  assumed  to 
be  present.  The  density  of  quartz  sand  is  assumed  to  be  1 .6  g/cc  which  corresponds  to  a  porosity  of 
27.4%,  assuming  an  initial  density  of  2.2  g/cc  for  solid  silica. 

For  dry  sand,  a  new  Equation-of-State  (EOS  7010)  for  fused  quartz,  which  was  recently 
developed  by  G.I.  Kerley  [14,  15]  was  used  in  this  study.  This  new  tabular  EOS  for  the  dry  sand 
pressure-density  relationship  gave  reasonably  good  agreement  with  Hugoniot  experimental  data  [16] 
for  dry  sand.  It  gives  correct  density  at  low  pressures  and  reproduces  the  phase  transition  in  the 
pressure  range.  It  is  clearly  a  better  model  than  the  EOS  7100  in  the  “seslan”  database  file  currently 
available  in  the  CTH  library  which  does  not  correlate  well  with  the  experimental  Hugoniot  data  for 
dry  sand.  Time  step  instability  problem  in  EOS  7100  was  drastically  reduced  in  the  new  EOS  model 
by  including  Maxwell  constructions  that  are  necessary  to  model  vaporization  behavior. 

The  porous  sand  is  not  expected  to  support  much  elastic  deformation  or  tension.  Therefore 
only  nominal  values  for  the  material  strength  and  fracture  parameters  with  a  yield  strength  of  1 .0E9 
dynes/sq.  cm,  Poisson’s  ratio  of  0.3  and  a  fracture  strength  of-1 .0E6  dynes/sq.cm. 

(b)  Wet  Sand 

For  a  porosity  of  27.4%,  fully-satiu-ated  wet  sand  would  have  a  composition  of  85.5%  silica 
and  14.5%  water  by  weight  and  a  density  of  1.87  g/cc.  An  EOS  table  for  this  composition  was 
constructed  by  G.I.  Kerley  [  1 4]  using  the  PANDA  mixture  model  [15].  The  procedme  was  similar  to 
that  used  for  dry  sand  except  for  the  addition  of  water  and  hydrogen. 

While  the  Hugoniots  for  wet  sand  (EOS  table  #701 1)  and  dry  sand  (EOS  table  #7010)  are  very 
similar  to  each  other  in  the  pressure-density  plane,  the  energy  absorbed  by  wet  sand  is  considerably 
lower  than  the  energy  of  dry  sand  at  the  same  shock  pressure.  Dry  sand  due  to  high  porosity  content  is 
able  to  absorb  more  of  the  explosive  energy  than  wet  sand  in  a  buried  mine  explosion.  As  a  result, 
An  explosion  in  wet  sand  more  acceleration  of  the  steel  plate  than  in  dry  sand. 

Due  to  high  water  content,  wet  sand  would  not  be  expected  to  support  significant  elastic 
deformation  or  tension.  Therefore,  the  yield  strength  of  wet  sand  was  set  to  zero.  Also  wet  sand  was 
not  allowed  to  fracture,  since  the  EOS  of  fully  saturated  sand  does  not  have  a  tension  region,  except  at 
very  low  temperatures  beyond  the  scope  of  the  cxurent  study. 
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RESPONSE  COMPUTATION 


The  simulation  of  the  propagation  of  the  explosive  blast  and  soil  products  through  the  air  and 
of  the  subsequent  interaction  with  the  steel  plate  was  performed  using  the  CTH  hydrodynamic  code. 
An  initial  time  step  of  1 .0  ps  was  selected  at  the  start  of  computation,  and  the  response  was  tracked  up 
to  2  ms,  corresponding  to  approximately  2,000  computational  cycles.  A  minimum  time  step  of  0.01 
ps  was  specified  to  avoid  numerical  instability. 

Gravity  effects  were  included  in  the  simulation  to  duplicate  actual  experimental  conditions. 
Active  Lagrangian  tracer  particles  were  selected  along  the  back  surface  and  the  horizontal  midplane  as 
well  as  the  vertical  axis  of  symmetry  of  the  plate.  The  selection  of  tracer  particles  where  history  plots 
of  position,  velocity,  pressure,  and  maximum  principal  stress  as  a  function  of  time  would  be  generated 
was  guided  by  experimental  considerations.  The  locations  selected  generally  corresponded  with 
locations  where  gauges  were  to  be  mounted. 

History  plots  of  vertical  linear  momentum  trapped  in  the  RHA  steel  plate  were  generated  for 
later  comparison  with  momentum  measurements.  Cross-sectional  plots  of  blast-soil-stracture 
interaction  and  target  response  were  plotted  at  firequent  time  intervals.  These  graphical 
representations  led  to  a  better  understanding  of  the  loading  mechanisms. 


RESULTS  AND  DISCUSSION 

Calculations  were  carried  out  with  the  Feb.  99  version  of  CTH.  The  modified  dry  and  wet  sand 
EOS  routines  provided  by  G.I.  Kerley  [13]  were  employed.  The  explosive  was  modeled  as  a  disk  of 
TNT,  30  cm  in  diameter  and  8  cm  in  thickness,  with  no  case.  The  initial  explosive  density  was  1 .60 
g/cc  and  the  corresponding  explosive  weight  was  9050  gm  (approximately  20  lb).  The  top  of  the 
explosive  charge  was  15  cm  below  the  soil  surface.  The  plate  was  244  cm  (8  ft)  square  and  20  cm 
t(approximately  8  in)  thick.  The  RHA  steel  was  46  cm  (  1 8  in)  vertically  upward  from  the  soil  surface. 
Additionally,  the  center  of  the  plate  was  offset  61  cm  (24  in)  horizontally  firom  the  center  of  the 
explosive  charge  as  shown  in  Figure  1. 

Cross-sectional  plots  of  blast-soil-plate  interaction  revealed  elongation  of  steel  plate  at  bottom 
comers  in  spite  of  high  elastic  modulus  and  stiffiiess  properties  of  the  RHA  steel.  Comer  deformation 
was  found  to  be  an  artifact  of  the  numerical  code  computation,  occurring  at  or  near  the  interface 
where  a  mesh  containing  a  relatively  stiff  material  is  located  next  to  a  soft  material.  This  condition 
caused  uncontrolled  material  flow  and  leakage  through  comers  along  the  interface.  It  was  decided  to 
use  the  boundary  layer  interface  (BLINT)  capability  in  the  code  along  the  interface  between  the  hard 
RHA  steel  target  and  the  relatively  soft  air  medium  to  mitigate  this  problem.  However,  no  significant 
change  was  observed  in  the  resultant  time  history  or  cross-sectional  plots  due  to  use  of  the  BLINT 
capability  available  in  the  hydrocode. 

Computation  using  the  20-cm-thick  anvil  plate  exhibited  some  transverse  bending 
deformation,  a  typical  cross-sectional  view  of  which  is  shown  in  Figure  2.  In  addition  to  bending 
deformation,  some  angular  rotation  was  observed  due  to  asymmetric  loading  fi'om  explosive 
detonation  in  diy  sand. 
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Figure  2.  A  Sectional  View  of  the  Mine-Soil-Plate  Interaction  at  2  ms  for  the  20-cm-Thick  Plate 

3D  Asymin.  mine-sand-plate  analysis 


Figure  3.  Comparison  of  Linear  Vertical  Momenta  Trapped  in  the  Steel  Plate  Due  to  Explosive 
Detonation  in  Dry  vs.  Wet  Sand  for  the  3-D  Asymmetric  Mine-Soil-Stmcture  Interaction  Problem 
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A  typical  history  plot  of  integrated  linear  momenta  in  the  vertical  Y  direction  for  the  steel 
plate  is  shown  in  Figure  3  The  momentum  is  initially  0  for  the  first  0.5  ms  until  the  arrival  of  the 
blast  wave  at  the  bottom  interface  of  the  plate.  The  momentum  absorbed  by  the  plate  increases 
rapidly  for  the  next  0.1  ms  and  then  gradually  for  an  additional  2.5  ms.  Thereafter,  it  reaches  a 
plateau  before  eventually  decaying  due  to  air  drag  and  gravitational  effects.  Comparison  of  the  anvil 
plate  vertical  momenta  from  detonation  of  a  mine  buried  in  fully  saturated  sand  (shown  as  a  curve 
joining  filled  diamonds)  vs.  dry  sand  (described  by  a  curve  joining  square  boxes)  indicates 
considerably  higher  momenta  for  wet  sand  relative  to  dry  sand  attaining  a  peak  differential  of 
approximately  150%.  This  is  significantly  higher  than  results  from  2-D  plane  strain  model 
computation  reported  earlier  [17  ].  It  appears  that  2-D  models  are  relatively  easy  to  set  up  and 
economical  to  run  and  the  results  are  quite  satisfactory  for  the  axisymmetric  case.  However,  for  the 
asymmetric  mine-soil-plate  interaction  problem,  simplified  2-D  models  are  inadequate  to  capture  3-D 
effects  necessitating  the  use  of  3-D  modeling  and  analysis. 

In  the  calculations  for  wet  sand,  the  soil  layer  does  not  thin  down  as  it  does  for  dry  sand.  This 
is  probably  because  wet  sand  is  not  porous  and  does  not  compress  as  much  as  dry  sand  when  shocked. 
Porosity  increases  the  Hugoniot  energy  at  a  given  shock  pressure,  so  that  porous  dry  sand  will  absorb 
more  of  the  explosive  energy.  Reduced  porosity  in  wet  sand  results  in  less  energy  absorption  and 
increased  contribution  to  the  impulse  delivered  to  the  plate.  However,  the  large  amount  of  momentum 
trapped  in  the  anvil  plate  for  the  wet  sand,  approximately  2.5  times  that  for  the  dry  sand  could  also  be 
contributed  by  the  sand  layer  in  front  of  the  expanding  explosive  reaction  products  which  is  thicker 
for  wet  sand  than  for  the  porous  dry  sand,  which  compresses  much  more  upon  shocking.  Also  wet 
sand  remains  cohesive  for  a  considerable  period  of  time  before  fragmentation  and  dispersion  occur  in 
the  debris  cloud,  while  lighter  sand  particles  have  low  mass  and  momentum  and  tend  to  separate  at 
early  times  Additionally,  the  increase  in  momentum  could  be  enhanced  by  increased  mass  due  to 
water  in  the  sand  and  increased  focussing  effects  of  wet  sand,  as  well  as  higher  stiffhess  and  shock 
speeds  in  the  wet  media.  Fmther  studies  will  be  needed  to  determine  which  mechanisms  have  the 
most  influence  on  structural  loading  due  to  explosive  detonation  in  wet  soil. 
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